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Abstract 

We describe a four-parameter model for non-exotic meson-meson scattering, which ac- 
commodates all non-exotic mesons, hence also the light scalar mesons, as resonances and 
bound states characterised by complex singularities of the scattering amplitude as a func- 
tion of the total invariant mass. The majority of the full S'-matrix mesonic poles stem from 
an underlying confinement spectrum. However, the light scalar mesons -ftrQ(830), ao(980), 
/o(400-1200), and /o(980) do not, but instead originate in ^Pq -barrier semi-bound states. 

In the case of bound states, wave functions can be determined. For cc and bb, radiative 
transitions have been calculated. Here we compare the results to the data. 
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1 Introduction 



Strong interactions have been under study for over half a century by now. Starting from curios- 
ity over nuclear forces |l|], they have been developing towards the need to understand Quantum 
Chromodynamics (QCD) 0. The modelling and analysis of meson-meson scattering is thereby 
believed to give valuable contributions. Here we will report on a model for meson-meson scat- 
tering which has been very successful in predicting and reproducing a host of experimental data 
with a very limited number of parameters 0-0. It is based on our perspective of QCD, as 
outlined below. 

Mesons are composed of quarks, antiquarks, and gluons, whether stable with respect to OZI- 
allowed hadronic decay, i.e., bound states like K, J/ip, and T, or unstable, decaying into two or 
more lighter mesons, like e.g. the resonances p, tu, and (f). Mesons can be described by systems 
of one constituent quark (g) confined to one constituent antiquark (g). Constituent quarks and 
antiquarks, which can be thought of as lumps of bare quarks surrounded by glue and a cloud 
of virtual qq pairs, determine the flavour of the meson and the bulk of its mass. The remaining 
strong interactions may be parametrised by an effective confining potential, e.g. a harmonic 
oscillator. 

Decay of a meson into a pair of lighter mesons is assumed to be the result of the creation of a 
new constituent qq pair with the quantum numbers of the vacuum i^Po), leading to OZI-allowed 
transitions. Mesonic resonances in meson-meson scattering are supposed to be formed through 
the inverse phenomenon, i.e., the annihilation of a constituent qq pair with the quantum numbers 
of the vacuum, one from each of the two mesons involved in the scattering process. 

The spatial quantum numbers (z/, J^*^ ) of a meson, either bound state or resonance, follow 
from the total spin of the constituent qq pair, their relative orbital angular momentum and 
their radial quantum number u. 
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Figure 1: Schematic picture of the model for the description of non-exotic meson-meson scatter- 
ing. 



In Fig. ([l|) we schematically show how the model for the description of non-exotic meson- 



meson scattering works. The employed dynamics is characterised by confinement and ^Pq qq-pa-ir 
creation. Confinement is provided by a harmonic oscillator potential with universal frequency uj. 
On the other hand, ^Pq gg-pair creation is modelled by a transition potential Vt, which has two 
parameters, that is, the universal overall transition intensity A, and a distance parameter a that 
is related to the shape of the transition potential. The so-called colour splitting of pseudo-scalar 
and vector mesons is absorbed in the transition potential through the parameter Qc, at least in 
the present version of the model. 

The four model parameters and the four constituent quark masses for non-strange (either u 
or d), strange, charm, and bottom quarks are adjusted to the data, i.e., to partial- wave cross 
sections of meson-meson scattering, or to mesonic bound-state and resonance spectra. Model 
predictions can then be obtained by giving as input data the relevant quantum numbers, namely 
flavour, isospin, total angular momentum, parity, and possibly charge conjugation. The outputs 
of the model are not only predictions for scattering quantities, but also comprise bound-state 
spectra and corresponding wave functions for the mesonic systems under study. 

The philosophy of the model is simple: in the scattering region for a meson-meson scattering 
process we assume the coexistence of a variety of different systems or channels, that is, one 
or more confined qq pairs, and all combinations of meson-meson pairs that have the correct 
quantum numbers. For instance, when studying scattering in an isosinglet, we assume that there 
are definitely non-zero probabilities to find in the interaction region an nn pair, an ss pair, but 
also a TTTT pair, a KK pair, and so forth. The first type of channels differs from the second type 
in that qq pairs are permanently confined, whereas meson-meson pairs feel more modest forces, 
which might even be neglected in lowest order. 

When we describe quarkonia by wave functions ipc and two-meson systems by wave functions 
il)f, then we obtain for their time evolution the wave equations 

(E-H,) ^,{r) = Vtijfif) and {E - Hf) (f) = [Vtf tfj^if) . (1) 

Here, He describes the dynamics of confinement in the interaction region, Hf the dynamics of 
the scattered particles, and Vt the communication between the two different sectors. 

For the dynamics of confinement we understand here that, as a function of the interquark 
distance r, the resulting quark-antiquark binding forces grow rapidly outside the interaction 
region. Consequently, we must eliminate ipc from the equations (|l|), since it is vanishing at large 
distances and thus unobservable. Formally, it is easy to do so. We thus obtain the relation 

(E - Hf) i,j (f ) = \Vtf (E - H,)-' Vt (r ) . (2) 

The set of coupled equations (|^) for the two- meson channels represented hy ipf can be solved 
for the scattering amplitude, either numerically or analytically, once He, Hf, and Vt are specified. 
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2 Fixing the parameters 



The two parameters A and a of the transition potential Vt can be fixed, together with the charm 
and bottom masses, by the spectra of the cc and bb vector states. A comparison of model 
results with experiment is shown in Table (|1]). In the first column of the table we find the 
spectroscopic notation for the various states: the first symbol indicates the quarkonium radial 
excitation N = u + 1, which for the here employed harmonic-oscillator confinement is related to 
the total bare quarkonium mass Eg^, by 

Eii, = uj (^u + ^ + 2) ^ quark masses . (3) 
In the last column of Table (|I|) we find the meson-meson channels which are involved in the 
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Table 1: Comparison of the model results for J^^ = 1 charmonium and bottomonium 

bound-state masses and resonance central mass positions to experiment p5|] . 



determination of the quarkonia spectra. Vector-vector channels come in pairs, one for total spin 
zero and one for total spin equal to two. Moreover, spin 2 may combine with F waves. The 
complete wave function of charmonium (bottomonium) is thus composed of twelve channels, 
the two cc {bb) channels, one for 5* wave and one for D wave, coupled to ten meson-meson 
channels through ^Pq transitions. Note that not all of these channels are open for the systems 
under consideration, but that does not necessarily mean their infiuence is negligible. The relative 
coupling intensities for the various channels, which follow from the three-meson vertices defined 
in Refs. [0, |ll|, are tabulated in Ref. |12|. In Nature, many more meson- meson channels couple 



in principle, but in the present version of the model only pseudoscalar and vector mesons are 
considered in the initial and final scattering states, since these are generally the ones that lie 
close enough to have an appreciable effect. 

The non-strange and strange constituent quark masses are fixed by the elastic P-wave vrvr and 
Kn scattering phase shifts, respectively, as shown in Fig. (^. We observe from the first plot that 
the agreement between model and experiment is not perfect for larger values of the vrvr invariant 
mass. We believe this to be a consequence of neglecting residual strong forces in meson-meson 
scattering that are not due to the ^Pq mechanism, like meson exchange. This might also explain 
why for Kit scattering the model gives superior results. Nevertheless, it is remarkable that a 
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Figure 2: Comparison of the model results (solid lines |^) for the tttt and Kn elastic scattering 
phase shifts with = 1~ quantum numbers to experiment (•), in the energy regions of the 
p(770) and i^r*(892) mesons, respectively. The data for tttt are taken from Ref. 0], and for Kn 
from Ref. Ma 



model which describes the cc and bb quarkonia bound states and resonances very well, can predict 
to a fair precision the tttt and Kir scattering data, just by fixing reasonable values for m„ and 
nis- Furthermore, we notice that the correct slopes of the phase shifts at 90°, as a function of 
the invariant masses, imply the model predicts the right widths for the p(770) and the i^'*(892). 

Finally, the parameter Qc, which accounts for the hyperfine mass splitting of pseudoscalar 
and vector mesons, is fixed by the mass differences between the irrjrj'K and pujipK* nonets. In 
Table @) we summarise the model results p for the light vector-meson sector. We find good 
agreement between model and experiment. However, with respect to the radial excitations of 



the p(770) meson, we believe that the particles refered to in Ref. [|15|, are the D states. In the 
harmonic-oscillator spectrum a D state is degenerate with the S state belonging to the next 
radial excitation. Through the coupling to the meson-meson channels the degeneracy is lifted, 
which gives rise to separate resonances (see Fig. |^ below). Also the 0(1680) corresponds in our 
model to the 0(1-D) state. 

For the pseudoscalars and the open-charm states, not shown here, the model predictions 
P] manifest much more disagreement with the data. Possible causes are a deficient treatment 
of chiral symmetry, relativity, and pseudothresholds. Further study on these points is clearly 
required, which is already in course. 
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state 


model exp. Ref. 


meson-meson channels involved 


p{lS) 
p{2S) 
p{lD) 
p{3S) 
p{2D) 


GeV GeV 
0.76 0.77 PI 
1.29 1.29 |17| 
1.40 1.42 |16| 
1.59 1.60 |18| 
1.68 1.68 |16| 


TTTT, KK, 
TTUJ, r]p, KK*, 
pp, K*K* 


(L,5) 

(1,0) 
(1,1) 

(1,0), (1,2) and (3,2) 


K*{IS) 
K*{2S) 
K*{3S) 


0.93 0.89 |15 
1.41 1.41 |15| 
1.73 1.72 |15| 


ttK, r]K, r]'K, 

TiK\ riK\ 7]'K*, Kp, Kuj, K(j), 
pK\ ujK% (PK* 


(1,0) 
(1,1) 

(1,0), (1,2) and (3,2) 


0(15) 
0(25) 
0(35) 


1.03 1.02 15 

1.53 

1.87 


KK, 

V0, KK*, 
K*K* 


(1,0) 
(1,1) 

(1,0), (1,2) and (3,2) 



Table 2: Real parts of the singularities in the meson-meson scattering matrices for some of the 
well-known light vector resonances. 
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3 Resonances 



The overall transition parameter A provides the communication between permanently closed qq 
channels and meson-meson scattering channels. To lowest order, this describes the decay of a 
meson, but to higher order also the occurrence of resonances in meson-meson scattering. For A = 
0, the set of equations (|l]) describes free, non-interacting two-meson systems with a continuum 
spectrum above threshold, and confined qq systems with a discrete mass spectrum, the so-called 
bare mesons, as given by formula (|^) for non-negative integer values of the radial quantum 
number u. When A is non-zero, then, for energies above threshold, the model describes the 
scattering of interacting two-meson systems, elastic as well as non-elastic, whereas, for energies 
below threshold, the model describes physical mesons that are stable with respect to OZI decays. 
A resonance does not represent just one single state, but rather a continuum of states. Each state 
under the resonance has a different mass, whereby the intensities of its couplings to the different 
two-meson systems are proportional to the corresponding partial-wave scattering cross sections 
at that energy. Hence, a resonance should be represented by a whole range of masses, though 
with varying importance for the scattering data. Such a phenomenon can also be described by a 
single complex number, namely a singularity in the analytically continued scattering amplitude. 
For narrow resonances, which here correspond to small values of A, the real and imaginary parts 
roughly equal the central resonance mass and half its width, respectively. However, for broad 
and highly deformed resonances, corresponding to the realm of strong interactions in our model, 
such simple relations do not exist. It would be helpful if in experiment complex masses could be 
measured, but exciting the imaginary parts of masses at particle accelerators goes far beyond our 
imagination as yet. At best, experiment supplies us with partial- wave cross sections, for which 
scattering amplitudes have to be constructed. Breit-Wigner shapes are the usual techniques, 
though for strong interactions probably not the most adequate tools. Moreover, since each 
scattering amplitude has a different pole structure, no generally applicable method exists how to 
locate singularities. Accordingly, since these poles, with the respective residues, constitute the 
most relevant information on scattering properties, it would be helpful if the Particle Data Group 
Collaboration could show us the world averages of full shapes of elastic and inelastic partial- wave 
cross sections or phase shifts, rather then just the central masses and widths of resonances. 

For the time being, it seems that we must content ourselves with at least agreement on the 
number of scattering-amplitude singularities within large intervals of the invariant mass. For 
small values of A, complex poles in the relevant Riemann sheet are expected to lie close to the 
bare meson masses, whereas for increasing A the singularities move farther away. Hence, one 
might expect that their number agrees with the number of bare mesons. Consequently, bound 
states and resonances could then both be labelled by the spectroscopic quantum numbers, as for 
instance shown in Table (HD. Nevertheless, in the case of the scalar mesons we will see below that 
not all singularities stem from the bare spectrum. 

Bare masses represent bare mesons that result from permanent confinement only. Therefore, 
models restricted to this particular sector of strong interactions should not adjust their param- 
eters to the bound states and resonances from experiment, although this is common practice 
unfortunately. In Fig. @ we depict how the real parts of the poles for the vector charmonium 
states behave in our model, when A is varied. For A = 0, the bare spectrum of formula (|^) is 
shown. For A = 1, one obtains the model result of Table (|^). Note that the effect of open- 
charm channels on the ground state is highly non-perturbative, despite of these channels being 
all virtual. 

Bound states can be represented by special singularities of the scattering amplitude, i.e., poles 
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Figure 3: The model results ^ for the central resonance positions of charmonium ^Si and 
^Di states, for four values of the parameter A, compared to the experimental situation |]15|. The 
various dashed lines indicate the threshold positions of the strong decay channels DD^ DD*, 
DsDs, D*D\ DsD:, and D^Dl 



on the real energy axis below the lowest threshold of the system. The J/tp (IS) is an example 
of such a state (see Fig. ^). It consists of two cc channels, a dominant one in 5* wave, the other 
one in D wave, as well as of channels with virtual pairs of D, Dg, D*, and D* mesons in P or 
F waves. The ip{2S) is an example of a bound state which originates in a bare meson located 
in the open-charm continuum. For values of A up to 0.6, the ip{2S) singularity of the scattering 
amplitude has an imaginary part, which corresponds to a resonance. But at the DD threshold, 
the ipi^S) pole turns real, for A ~ 0.6. For the model value of A = 1, the singularity is found on 
the real axis below the DD threshold, in accordance with the physical V'(2S) state. 
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4 Scalar mesons 



In the foregoing, we have discussed how the parameter A, which describes the overall intensity 
of the coupling of qq systems to meson-meson scattering channels through the ^Pq mechanism, 
can be adjusted to the data. When we carefully inspect wave equation (^, then we observe that 
the centrifugal barrier is the main contribution of the left-hand side of the equation, whereas 
the right-hand side is proportional to A^. Consequently, for S-wave scattering, in the absence 
of a centrifugal barrier, the right-hand side of Eq. (|^) dominates the interaction. Now, for low 
energies far away from the bare states, this term has the form of a potential barrier, in which 
loosely bound two-meson systems may exist, solely depending on the value of A. This is exactly 
what happens for the light flavours. The ^Pq barrier, not high enough to form bound states, gives 
rise to a complete nonet of light scalar resonances with central mass positions at about 0.8-1.0 
GeV. 

bare scalar first radial second radial 

qq nonets ground state excitation excitation 

• • • 

« lAGeV ^ l.SGeV « 2.1GeV 



S-wave meson-meson V 
resonance nonets 



coupling 
to meson loops 



V 



V 



O.SGeV 
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Figure 4: Schematic presentation of pole doubling for the light flavours [^]. Bare meson nonets 
turn into nonets of mesonic resonances. In the process of switching on the model parameter A, one 
extra nonet of singularities appears in the nonet of S'-wave meson-meson scattering amplitudes. 

For S'-wave meson-meson scattering, our model produces, besides an infinity of = 0^ 
5'-matrix poles or resonances, each stemming from a level in the bare spectrum, also one ad- 
ditional singularity that is of a completely different nature. In Fig. (^), we schematically de- 
pict this situation. The nonets of scalar resonances are formed by P-wave qq states coupled 
to pseudoscalar-pseudoscalar and vector-vector two-meson channels. There are two isodoublets 
(non-strange-|-strange), which among their decay modes have a Kit channel. Then there is one 
isotriplet (non-strange+non-strange) involving an r^vr channel. And finally there is a system of 
two isosinglet channels, one for nn and one for ss, which mostly couple to each other through 
their common tttt and KK channels. 

Each of the isodoublet states has an extra pole at 727 — i263 MeV, which corresponds 
to the two doublets of i^Q(830) resonances. The three isotriplet states have extra poles at 
962 — i28 MeV, which describe the triplet of ao(980) resonances. Finally, the coupled system of 
isosinglets has extra poles § at 470-i208 MeV and 994 -z20 MeV, which yield the /o (400-1 200) 
and /o(980) resonances, respectively. 
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The wave functions of the K*{830), ao(980), /o (400-1200), and /o(980) resonances do not 
have nodes. Consequently, in this sense they must be considered ground states. The ao(980) 
and /o(980) resonances are experimentally well established. The /o(400-1200) resonance is still 



disputed, whereas the i^'o(830) has recently been reported in preliminary experimental work pO |. 

The nonet of resonances stemming from the nonet of ground states of the bare meson spectrum 
around 1.4 GeV consists of the experimentally estabhshed scalars -^"0(1430), ao(1450), /o(1370), 
and /o(1500). The higher nonets are still incomplete. 

Scalar mesons, especially the light ones, have been the subject of lively debate for several 
decades. We believe that the here described model has all the necessary ingredients for a rea- 
sonable description of the data (see also Ref. PT|). 
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Figure 5: Elastic S-wave Kn cross section. Data are taken from Refs. [Tj, (0) and p3| 



The solid line is a model result Maximum of the theoretical cross section is at 832 MeV. 

In Fig. d^), the prediction of Ref. [|] for the 5"- wave cross section in isodoublet tiK scattering is 
shown. One observes that the central peak lies at 832 MeV. Nevertheless, the related singularity 
of the irK scattering amplitude we found at 714 — ^228 MeV. Hence, straightforward relations 
between bumps in the cross section and pole positions do not always exist. 

Recently it has been recognised that S'-matrix pole positions have more consistency than 
resonant structures in the data. Values found in the literature for the a and k, resonances are 
collected and depicted in Table (0). The world averages for the /o (400-1200) and Kq{830) pole 
positions, thereby ignoring the differences in the nature of the various analyses, are 

/o [(506±119)-i(246±85) MeV] ; K* [(786 ± 78) - i(255 ± 96) MeV] . 
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date Ref. 



a (MeV) 
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1973 
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910 


— 


i350 






420 


— 


i370 






370 
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506 


— 
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387 


— 


z305 






470 
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i250 
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559 


— 


il85 






469 


— 


il79 






602 


— 


il96 


875 - 


i335 


442 


— 


z227 


770 - 


i250 


422 


— 


z213 






445 


— 


z235 






420 


— 


i317 






518 


— 


i261 






445 




i221 


779 - 


i330 








911 - 


il58 


498 




i321 






732 




il23 






478 




il62 


797- 


i205 








708 - 


i305 


449 




i241 






585 




il93 


905 - 


i273 


535 




^55 












714- 


i228 



K (MeV) 



0.4 0.6 0.8 GeV 



-0.1- 


• 




• 


-0.2- 
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• 

••• 


• 




• 






-0.3- 
-0.4- 


• 

• • 
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-0.5- 
GeV 


• 

Im[^/s) 




CT 










0.7 

1 




0.8 GeV 

1 , 1 








• 


-0.1- 
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-0.2- 


• 


■ 


• 


-0.3- 




• 


• 


• 

• 


-0.4- 


• 






GeV 


Im{yfs) 




K 



Table 3: Results and references for a and k poles in the complex ^/s plane for isosinglet tttt 
and isodoublet t^K S-wave scattering, respectively. The upper figure shows the locations of the 
a poles, the lower one for the k poles. World averages are indicated in the figures. Our model 
predictions are indicated by o. 
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5 Radiative transitions 

Multichannel wave functions are neither very suitable for figures, nor for comprehensive tables. 
Radiative decays, however, represent a rather sensitive test for their correctness. Such tests has 
been carried out in Ref. for the cc and bb sectors of the model. It involves the development 
of a formalism for radiative transitions of multicomponent systems with quarkonium and meson- 
meson channels. In Tables (^) and (|^), we compare the theoretical predictions of Ref. for 
the cc and bb systems, respectively, with the presently available data. 



state 


decay product 


model |51 


experiment 15| 






keV 


keV 


jms) 


IVcilS) 


1.67 


1.13±0.41 




7^co(lP) 


31.7 


25.8±5.4 




7^ci(lP) 


49.8 


24.1±4.9 




7^c2(lP) 


38.5 


21.6±4.6 






1.41 


0.78±0.25 






52 


98±43 


X,i(lP) 




210 


240±52 






228 


270±46 



Table 4: Comparison of the model results for cc radiative transitions to experiment. 



state 


decay product 


model |51 




experiment 


m 






keV 


keV 


T{2S) 




1.6 


1.7±0.5 




7^w(lP) 


2.1 


3.0±0.8 






2.6 


3.1±0.8 


T(3^) 


7^w(2P) 


1.6 


1.4±0.3 




7^61 (2P) 


2.6 


3.0±0.6 






3.1 


3.0±0.6 



Table 5: Comparison of the model results for bb radiative transitions to experiment. 

The agreement of the model's predictions with the data is fairly good and has much improved 
since their publication in 1991, which is an additional argument in favour of the usefulness of the 
employed formalism. 
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6 Data analysis 



The here described model produces non-exotic bare mesons with all possible angular quantum 
numbers and an infinity of radial excitations, which, through the universal parameters A and a, 
turn into the experimental spectra of all existing mesonic bound states and resonances. However, 
in spite of the model's success to reproduce a host of experimental data with a very limited number 
of parameters, it is clearly not suited as a tool for data analysis, owing to the specific model choice 
of the confining qq potential, and furthermore the rather complicated matrix expressions needed 
to obtain S-matrix-related observables. 

In recent work ||^, we have indicated how the model's formalism may be applied to data 
analysis. The bare meson masses and their couplings to the meson-meson channels enter as free 
fit parameters in this approach. In particular, we have shown that the ^-matrix singularities 
describing the light scalar mesons also appear here. Hence, this simplified model, which can be 
solved analytically, serves well as a laboratory for studying the general properties of meson-meson 
scattering. Moreover, for a more detailed description of specific processes, parts of the analytic 
solutions can be substituted by adjustable parameters. 



7 Hybrids 



In a recent work |^|, Frank Close suggests the use of a model similar to ours for the study of 



hybrids. Certainly, as the objects which we here call constituent quarks can exchange quark and 
colour degrees of freedom, one may naturally think of different quark and/or colour configura- 
tions. In this sense, also exotic mesonic resonances could, in principle, be accommodated. 
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8 Conclusions 



Bound states and resonances are most conveniently characterised by singularities of the scattering 
matrix as a function of the total invariant mass. It has the advantage that even in cases where 
the resonating behaviour of the scattering cross sections is not very clear from experiment, 
their existence clarifies the classification of well-established resonances into flavour multiplets . 
However, a disadvantage is that each different theoretical model locates such poles at different 
complex values. Hence, no agreement can be expected on their positions. Nevertheless, when 
several analyses, describing the same phenomenon, obtain poles at more or less the same values, 
like those given in Table (^, then no doubt can exist about their existence in the true model of 
Nature. 

Knowledge about how the poles depend on the parameters of strong interactions in each 
model could, moreover, lessen the discrepancies. Since it is indispensable to gain experience on 
the behaviour of mesonic S'-matrix singularities, it is useful to dispose of analytical expressions 
that incorporate most of our knowledge on strong interactions, and are also capable of describing 
the data to a moderate degree of accuracy. The outlined simplified model serves well for that 
purpose. 
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